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Using the transformation theory for the Heun equation, we find substitutions of variables which 
preserve the form of the equation for the kernels of integral relations among solutions of the Heun 
equation. These transformations lead to new kernels for the Heun equation, given by single hyper- 
geometric functions (Lambe- Ward-type kernels) and by products of two hypergeometric functions 
(Erdelyi-type) . Such kernels, by a limiting process, also afford new kernels for the confluent Heun 
equation. 
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C"| ■ 1. Introductory remarks 



The group of transformations of variables which changes Heun's equation into another version 
of itself was initially established by the very Heun in 1889 [l[ and fully accounted in 2007 by 
Maier 0], who tabulated the 192 substitutions in detail by writing explicitly the transformations 
of each parameter and variable of the equation. Firstly, we show that it is possible to construct 
a similar table for the kernels of integral relations among solutions of the equation. 

In the second place, we show that some of these transformations generate new kernels given 
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^ c*j by hypergeometric functions when applied to the kernels found by Lambe and Ward in 1934 
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and new kernels in terms of products of hypergeometric functions when applied to the kernels 
found by Erdelyi in 1942 Finally, by means of a limiting procedure we get new kernels also 
for the confluent Heun equation (CHE) . 

The transformations of the Heun equation and its integral relations suppose the usual alge- 
■ braic form for the equation [l|, 0, 0] , namely, 
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x(x — l)(x — a) 



H = 0, [e = a + (3+l-j-5] (1) 



where a G C \ {0, 1} and x = 0, 1, a, oo are regular singular points with indicial exponents given 



CN ' by {0, 1 — 7}, {0, 1 — 5}, {0, 1 — e} and {a, /?}, respectively. The constants a, a, f3, 7 and <5 arc 



called singularity parameters, whereas q is called accessory parameter since it is not associated 
with the singular points or their indicial exponents. 
By defining the operator M x as 

d 2 d 
AI X = x(x — l)(x — a)-— it + [7(2; — l)(x — a) + 8x(x — a) + exix — 1)1 — — H a(3x (2) 
ax ox 

and, by interpreting this as an ordinary derivative operator, the equation reads 

[M x - q] H{x) = 0, [a^Oorl], (3) 

The invariance of the equation with respect to the replacement of a by /3 does not imply that 
its solutions are symmetric in a and /3; it simply means that the substitution of a for j3 leads 
to another solution. The values a = and a = 1 arc excluded because in these cases there are 
only three singular points and then the equation may be reduced to the Gauss hypergeometric 
equation 

d 2 F dF 

u(l-u)-—+ [c-fa + b+lWl- akF = 0, (4) 

du l du 

where u = 0, 1, 00 are regular singular points with indicial exponents {0, 1 — c}, {0, c — a — b} 
and {a, b}, respectively. 
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By keeping a, 7 and S fixed, the confluence procedure is given by the limits 

a, /3, q — > 00 such that — — > > —p, — — > —a, (5) 

a a a 

where p and a are constants. This yields the CHE or generalised spheroidal wave equation 
1, 

H dH 
x(x - l)-^2" + [ - 1 + (l + S)x + px(x - 1)] — + [apx - cr]H = 0, (6) 

where x = and a; = 1 are regular singularities, whereas x = 00 is an irregular singularity. 
On the other hand, "K(x) is defined by 0-0] 



M(x)=/ w(x,y)G(x,y)H(y)dy = y^-\l - y) 5 ~ x (l - y -) G(x,y)H(y)dy, (7) 
•Ah -A/i a/ 



where ff(a;) represents a solution of equation ([TJ. Then, J£(x) will be a solution of the Hcun 
equation if: (i) the kernel G(x, y) is solution of the partial differential equation 

[M x -M y ]G(x,y)=0, (8) 

where M y is obtained by setting x = y in the expression for M x , (ii) the integral exists and 
(iii) the limits of integration arc so chosen that the bilinear concomitant P(x,y), given by 



P(x,y)=yi(l-y) s l-£ 



dG(x,y) ,dH(y) 
H (y) — 5 G(a;,y)— — 



(9) 



fulfills the condition P(x, j/i) = P(a;,2/2)- In Appendix A we show how these equations are 
obtained from the general theory of integral relations @ ■ 

By the choice given in Eq. ([7J for the weight function w(x,y), equation ([5} for the kernels 
is expressed in terms of the operator M x which appears in the Heun equation ([T]) and in terms 
of the functionally identical operator M y obtained by setting x = y in M x . Then, in order to 
establish the transformations of the kernels it is sufficient to demand that M x and M y transform 
in the same way. For the Heun equation these transformations will be inferred from the Maier 
transformations for the Hcun equation. 

By using suitable weight functions the above result holds also for the other equations of the 
Hcun family, that is, for confluent, double-confluent, biconfluent and triconflucnt Hcun equa- 
tions. Then, the transformation for the kernels may be inferred from the known transformations 
of each equation ||. It seems that this connection has not been explored as yet [9MTTI| . 

However, the transformations become effective only if we know an initial kernel. For the Heun 
equation, new kernels in terms of single hypergeometric functions will be generated from the 
kernels found by Lambe and Ward Q , while kernels given by products of two hypergeometric 
functions will arise from the ones found by Erdelyi Q ■ These afford initial kernels for the CHE 
by the limiting process ((SJ). In addition to kernels given by confluent hypergeometric functions, 
we find kernels given by hypergeometric functions, products of two confluent hypergeomet- 
ric functions, and products of one confluent hypergeometric function and one hypergeometric 
function. 

In section 2, firstly we present the 8 so-called index or nomotopic transformations which do 
not change the independent variable x, and the 24 Mobius or homographic transformations 
which result from linear fractional substitutions of the independent variable. Composition 
of such substitutions gives the group of 192 transformations. After this, the transformations 
are extended to the kernels of the equation, and these are used to generalise the kernels of 
Lambe- Ward and Erdelyi. 

In order to generate the full group by composition of homotopic and homographic trans- 
formations, it is necessary to use the index transformations in Maier's form. This remark is 
important for avoiding incorrect results. For example, the forms given in Refs. Q and [T2I ] are 
inappropriate as we shall explain in section 2.1. 

The kernels for the CHE are obtained in section 3, where we introduce as well the transfor- 
mations of Eq. ^ and its kernels. In section 4 we point out that even for the double-confluent 
Hcun equation (DCHE) it is possible to determine new kernels by using again a limiting pro- 
cess, and discuss how to transform certain solutions of the Heun equation into solutions useful 
for applications. Appendix A provides a derivation of Eqs. ([8]) and ©, while Appendix B lists 
the Mobius transformations for kernels of the Hcun equation. 
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2. Heun's equation 



First we examine the transformations for the Heun equation, emphasising that it is not 
allowed to permute the parameters and ft in the homotopic (index) transformations. Second, 
we obtain a general prescription for the transformations which preserve the equation for the 
kernels and write, explicitly, the index transformations for the kernels. In the third and fourth 
subsections, respectively, we generalise the kernels given by hypergeometric functions and by 
products of hypergeometric functions. 



2.1. Transformations of Heun's equation 



There are 24 (including the identity) Mobius substitutions of the independent variable x 
which leave the form of Heun's equation invariant, in general after a change of the dependent 
variable. They are given by fractional linear transformations x n> g(x) = (Ax + B)/(Cx + D), 
AD ^ BC, which map three of the points 0, 1, a and oo onto 0, 1, oo. The expressions for g(x) 
are displayed in the matrix 
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(a—l)x 
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a(x — 1) 



^ „ x — 1 a — x . 1 x — 1 x — a 1 x — a a 

x — a a ' x x x 1—x x — 1 a—x 

1-x a(x-l) g—x . a a(x — l) x—a 1-a x—a a-1 

1 — a x — a a— 1 ' x (a— l)x (1— a)x 1—x a(x— 1) a—x 



(10) 



where the elements in each column are proportional to one another and, in each row, the 
elements after the semicolon are the inverses of the elements before semicolon. For the identity, 
and for (x — a)/(x — 1), a(x — l)/(x — a) and a/x the other singular point is mapped onto a, 
while for the remaining cases it changes to [U, [j| 

1 1 a a — 1 , . 

-, 1-a, z , r, ■ (11) 

a 1 — a a — 1 a 

Sometimes solutions for the Heun equation arc denoted by H(x) = Hl(a,q;a,0,'y,5;x), 
where HI means 'Heun- local', that is, a solution which converges in a region containing 
only one of the four singular points 0, @- F° r brevity, we drop the letter I, writing 
H(x) = H(a, q; a, /?, 7, S; x). Then, the Mobius substitutions permit onto map a solution H(x) 
into new solutions according to 

H(a,q;a,0,'y,S;x) \-> f(x)H a,q;a,j3,j,5;g(x) , (12) 

where the prcfactor f[x) symbolises the transformation, if any, of the dependent variable which 
brings the differential equation with the variable g(x) into a Heun equation having parameters 
a, q, a, (3, 7 and 6. Depending on the transformation considered, we have 

f{x) = l, x- a , (l~x)- a , [l-x/a)- a ,ov (13) 



f(x) = l, x-e, (i-x)-e, [l-x/a)-^ (14) 

up to a multiplicative constant. The prefactor f(x) = 1 corresponds to the linear transforma- 
tions, namely: g(x) = 1 — x, (a — x)/a, x/a, (1 — x)/(l — a) and (a — x)/(a — 1). The first form 
(fl3|) is the one that will be adopted in the present article. 

On the other side, the index transformations do not change the independent variable. They 
are given by 8 elementary power transformations of the dependent variable @, [B[ , namely, 

H(a, q; a, /3, 7, 5; x) H> x Tl (1 - x) T2 (1 - x/a) Ta H(a, q; a, (3, 7, 5; x), (15) 

where t\, ti and T3 are the indicial exponents at 0, 1 and a, respectively, namely: n = or 
1 — 7, T 2 = or 1 — S, and T3 = or 1 — e. Since there is no change of the independent variable, 
the positions of the singular points remain fixed, in contrast with the fractional transformations. 
For this reason, they are also called homotopic transformations. 

The composition of these two types of transformations (elementary powers and fractional) 
generates the group containing the 192 transformations given in Maier's table We refer to 
such transformations by Mj (i = 1, 2, • • • , 192) following the order in which they appear in the 
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table, Mi being the identity transformation. Thence, by regarding the M, as operators, the 
effects of both transformations on a solution H{x) are represented by 

MiH(x) = MiH(a,q;a,P,y,S;x) = fi(x)H [at,qi;ai, /?i,7i,&; , (16) 

where is given by a or one of the expressions written in ([lip , whereas Qi(x) is one of the 
elements of the matrix (fTOj) . 

In Maier's table the elementary power transformations appear in the entries Mi — M4 and 
M25 — M28, but here they are denoted by T (£ = !,••• ,8) according to the correspondence 



Ti T 2 T 3 T 4 T 5 r 6 T 7 T 8 
Mi M 25 M 2 M 26 M 3 M 27 M 4 M 28 



(17) 



where T\ is the identity: T±H(x)=H (a, q; a, /3, 7; <5; x). The transformations T 2 , T 3 and T5 are 
given by 

T 2 H{x) = x^H [a, «? - (7 - l)(<5a + e); ,3 - 7 + 1, a - 7 + 1, 2 - 7, 5; a] , 

r 3 #(a;) = (1 - x) 1 -^ [a, q-(S- 1)70; - 5 + 1, a - <5 + 1, 7, 2 - 6; x] , (18) 

T 5 H(x) = [1 - (or/a)] 1_e ff [a, g - 7(0* + /3 - 7 - 5); -a + 7 + J, + 7 + 5, 7, 5; x] . 

and are the generators of the other Ti. In these transformations we cannot change the order of 
the parameters cti and ft, that is, we must read 

for T 2 : a 2 = /? — 7+I, ft> = a — 7+1; 

forT 3 : a 3 = /?-<5+l, /? 3 = a-£+l; (19) 
for T 5 : a 5 = -a + 7 + 5, ft = -/3 + 7 + (5. 

In effect, it is possible to obtain the 192 transformations only if T 2 , T 3 and T5 transform the 
prcfactors x~ a , (1 — x)~ a and (1 — x/a)~ a of the Mobius transformations into x~@ , (1 — x) _ ^ 
and (1 — x/a)~P , and vice-versa. For the other transformations, the positions of a, and ft 
result from the compositions 

T i H{x)=T 2 T 3 H{x), T e H(x) = T 2 T 5 H, T 7 H(x) = T 3 T 5 H(x), T S H = T 2 T 3 T 5 H(x), (20) 

where the order of the operators T is irrelevant on the right-hand side since TTj=TjT. In 
spite of this, for the three transformations Sleeman and Kuznetsov [l2| writes c*i and Pi in the 
inverse order, while Arscott [f| inverts the order in T 2 and T 3 . 

The above order for a,i and ft is important regardless the composition among the two types 
of transformations. This becomes apparent by considering an Erdelyi solution in series of 
hypergeometric functions F(a, b;c;x), given by [l3[ 



00 

Hi(x) = V F{n + a,-n-a-l+j + S; T} x), (21) 



n=0 



where the coefficients fo^ 1 ' satisfy three-term recurrence relations. From this we can generate 
a subgroup constituted by 8 solutions by writing Hi(x) = TiH\{x). In particular, the solution 
with a in the place of /? is 

00 

H 3 (x)=T 3 H 1 (x) = (l-x) 1 - 5 ^^? F(n + l3 + l-S,-n-f3 + T}T ,x) 

n=0 

00 

= F(n + /?,-n-/?-l + 7 + <5;7;*), (22) 



where the last equality follows from Eq. (|32p written later on. By interchanging a 3 and ft in 
T 3 , we would obtain the identity H 3 (x) = H\{x) , that is, we would miss one solution at least. 

As aforementioned, we take the Mobius transformations with prefactors given in Eqs. (|13| 
as the basic ones. Then, the substitutions (fTU)) correspond to the following entries in Maier's 
table: 



Mi M 5 M13 M 49 M57 M101; M 145 M 53 M 97 M 149 Mi 05 M 157 
M 9 Mai Mi 7 M 6 i Af 65 Mi 09 ; Mi 53 M 69 M 117 M W1 M 113 M W5 



(23) 



where each column presents the same prefactor. 
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2.2. Transformations of kernel equation and notations 

Since Eq. ((Sj is independent of the parameter g, a general kernel will be denoted by G(x, y) 
= G(a; a, /3, 7, <5; x, y). Then, the transformations which take the place of previous Mobius and 
index transformations are given, respectively, by the mappings (symmetrical in x and y) 



G [a; a, /3, 7, 5; x, y] h-> f(x)f(y)G a; a, /3, 7, 5; p(x), £>(y) 
G [a; a, /?, 7, <5; x, y] ^> [xyp [(1 - x)(l - y)p 







T3 




[('-:) 


G 


a 











a; a, £,7, 8; x, y 



where the prefactors f(x) and f(y), as well as the fractional transformations q{x) and g(y), are 
formally the same which occur in the transformations of the Hcun equation. These substitutions 
preserve the form of equation [M x — M y ]G = for the kernels because all the parameters of the 
operators M x and M y transform as in the Heun equation, and constant terms corresponding to 
the transformations of q cancel out. In terms of operators we rewrite these transformations as 

3QG(x, y) = XiG(a;a,f3,-f,5;x,y) = fi(x)fi(y)G [<n; on, fa, 7,, Si\ Qi(x), Qi(y)\ , (24) 



where %i is obtained from the corresponding Mi of Eq. (pi 

We split the operators %i in two subgroups: the operators denoted by Ni correspond to 
the 8 index transformations Tj , and the operators denoted by Ki correspond to the 24 Mobius 
transformations. Thus, the nomotopic transformations Aj for the kernels are 

NxG(x,y) = G(x,y) = G [a; a, /3, 7, 5; x, y] [Identity], 

N 2 G(x,y) = (xy) 1 -^G[a;P- 1 + l,a- 1 + l,2- 1 ,S;x,y]. 

N 3 G(x,y) = [(l-ar)(l-y)] 1 - 4 G[o; J 8-5 + l,a-<y + l > 7,2-5;x > i/] > 

N 4 G(x,y) = (xyY-^l-xXl-y} 1 - 6 

x G [a; a - 7 - S + 2, /3 - 7 - 6 + 2, 2 - 7, 2 - 5; x, y] , 



N 5 G(x,y) 



l- X - 

a 



1- V - 



G [a; -a + 7 + 5, -(3 + 7 + 5, 7, (5; x, y] , 



N 6 G(x,y) = (xy) 



1-7 



DO 



G[a; -/3 + 5 + 1, -a + 8 + 1, 2 - 7, 6; x; y] 



N 7 G(x,y) = [(l-x^l-y)] 1 - 5 

x G[a;-/3 + 7 + l, -a + ~f + l,~/,2 - 6;x,y], 



1-*) (l- V - 
a) V a 



l-e 



A 8 G(x,y) = (^^[(l-xXl-y)] 1 - 5 

x G[a;2-a,2-/3,2-7,2-<5;a;,j/] 



l- X -\ (l-V. 
a) V a 



Notice that in fact Ai is the identity only if a 1— >■ a and /3 m-/?. As in X^, we cannot 
permute a* and fa in A2 , A3 and A5 because such transformations must change the exponents 
a of the prefactors (|27p (for fractional transformations) into fa and vice-versa. For the other 
transformations, the above positions of cti and fa result from the compositions 



N 4 G(x,y) = N 2 N 3 G(x,y), 
N 7 G(x,y)=N 3 N 5 G(x,y), 



N 6 G(x,y) = N 2 N 5 G(x,y), 
N s G(x,y) = N 2 N 3 N 5 G(x,y), 



(25) 



where the order of the operators Ni is irrelevant on the right-hand side. Thus, A2, A3 and A5 
are the generators of the index transformations for the kernels. 
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On the other side, each of the 24 Mobius transformations Mj given in the matrix (|23l) is 
associated to a kernel transformation denoted by Kj (j = 1, • • • ,12 for the first row and 
j = 13, • • ■ , 24 for the second) as 



K x K 2 K 3 K 4 K 5 K 6 ; K 7 K 8 K 9 K w K u K 12 
K13 K15 Kia Kn Kig; K19 K20 K21 K22 K23 K24 



(26) 



As in the homotopic transformations, first we must find the transformation Mj for H(x) in 
Maier's table and, then, write the kernel transformation Ki by using the prescription given 
in Eq. (|24[) . The 24 expressions for Ki are written down in Appendix B and will be used in 
following subsections. The prefactors for those transformations are 

1, (zy)- a , [(l- X )(l-y)}-\ [(l-£) (l-£)p- (27) 

Transformations having prefactors with exponent j3 are generated from the ones of Appendix 
B by applying N 2 when the prefactor is (zy)~ a , A 3 when [(1 — x)(l — y)]~ a and A 5 when 
[(l-x/a)(l~y/aT a . 

We could denote an initial kernel by G 1 (x,y), and by G^(x,y) the kernels obtained from 
G\(x,y) by using Ni and Kj, where 



j indicates Mobius transformation, Kj, 
i indicates index transformation, TV,-, 



being necessary to specify the transformation applied in the first place since in general index 
and Mobius transformations do not commute. However, this notation is not sufficient because 
we will consider a set of initial kernels rather than a single kernel. For the Lambe-Ward case, 
we start with six kernels given by distinct hypergeometric functions and use the notation 
where k runs from 1 to 6; for the Erdelyi case we take 36 products of hypergeometric functions 

l(k I) 

and then the initial set is denoted by G x ' . Thence, the actual notation will be 

Gj (x,y) for Lambe- Ward-type kernels ; G^' 1 ' {x,y) for Erdclyi-type, (28) 

where the indices inside parentheses are not affected by the application of the transformations 
iV, and Kj. 

The Lambe-Ward kernels G^ k \x,y) defined in Eq. (|34|) and Erdelyi kernels G]^ k ' l \x,y) 
defined in Eq. (|47[) are the initial kernels which are obtained by solving directly the kernel 
equation. We will find that Kj with j = 2, ■ • • ,6 are the only effective Mobius transformations. 
These lead to five additional sets of Lambe- Ward-type kernels, denoted and obtained as 



G 2( fe ) = K2G l(k) G 3« = KaG l (fc) G 4( fe ) = ^(k) 

Gf k) =K b G\ (k \ Gf k) =K 6 G\ {k \ [fc = l,2,-",6]. 
Similarly, the new Erdclyi-type kernels arc 

r-2(k,l) r , fl(k.l) /-3(fc,Z) T r r-Uk.l) r-Uk.l) r , fl(k.l) 

G? (M) =K 5 G{ ik ' l) , Gl™ = K a G\M , [k,l= 1,2,- •■ ,6]. 



(29) 



(30) 



In each case, the subscript could assume eight values when we apply the homotopic transforma- 
tions Ni. Nevertheless, we will find that only three of the TV, are effective due to fact that one 
of the generators A^, A3, A5 becomes equivalent to the identity or two of them are equivalent 
to each other. Thus, there are only four values for the subscripts. 



2.3. Generalisation of Lambe- Ward's kernels 

The Lambe-Ward as well as the Erdelyi kernels are given by hypergeometric functions 
-F(a, b;c;u) = F(b, a;c;it). In fact, in the vicinity of the singular points 0, 1 and 00, the 
formal solutions for the hypergeometric equation (|4]) are Q , respectively, 



FW(u)=i?(a,b;c;u), F^{u) = u^F (a + 1 - c, b + 1 - c; 2 - c; it) ; (31a) 
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^(3) ( u ) = F (a, b; a + b + 1 - c; 1 - it) , 

FW (u) = (1 - u) c ~ a ~ b F (c - a, c - b; 1 + c - a - b; 1 - u) ; 



(31b) 



(31c) 



F^iu) = u" a F (a, a + 1 - c; a + 1 - b; 1/u) , 
F( 6 \u) = u~ h F (b + 1 - c, b; b + 1 - a; 1/u) . 
Each of these functions may be written in four forms by using the relations 

F(a, b; c; u) = (1 - u) c - a - b F(c - a, c - b; c; u), (32) 

F(a, b; c; u) = (1 - u)~ a F [a, c - b; c; u/(u - 1)] . (33) 

The Lambe and Ward kernels Q in terms of single hypergeometric functions are obtained 
by taking u = xy/a and G(x, y) = F(u) 

in [M x — M y ] G(x, y) = 0. Hence, F(u) satisfies equation (0| with a = a, b = /? and = 7. 
Thus, the six formal kernels have the form 

G} (I) (x, y) = F« (^) , [a = a, b = ft c = 7] (34) 

In this manner, up to a multiplicative constant, the initial set of kernels is given by 



G^'(x,y) = F[a,P;r,- 

a ' (35a) 
G} (2) (x, y) = {xyf-^F (a + 1 - 7, + 1 - 7; 2 - 7; ^) ; 



Gj (3) (a, y) = F (a, ft a + /3 + 1 - 7; 1 - ^) , 

(x, y)=(l- f ) ^ F (7 - a, 7 - ft 1 + 7 - a - ft 1 - f ) 



Gj (x, 1/) = {xy)- a F a,a + l- 7 ;a + l- /3 
Gi (6) (x, y) = (^)-^F (p + 1 - 7, ft P + 1 - a; — 



(35c) 



By using this set of initial kernels, some of the kernel transformations become superfluous. In 
effect, by the transformations iVj we could obtain a subgroup containing eight sets. However, 
for the present case N2 is ineffective since 



N, 



/ r l(l) r l(2) -1(3) r l(4) -1(5) r l(6)\ _ ( r l(2) -1(1) -1(3) r l(4) -1(6) r l 



(5) 



that is, N2 simply rearranges in a different order the previous kernels. In this sense the generator 
N2 is equivalent to the identity N\ and, so, the index transformations can generate only three 
additional sets due to composition relations (f25|) . In fact we find that 



Therefore, it is sufficient to use the transformations N3, N$ and N7 to produce three additional 
sets, namely: G3 , Gg and G^ k \ The eight kernels G* and G\^ (i = 1,3,5,7) coincide 
with the ones given by Lambe and Ward. Now, we regard the generalisations arising from the 
Mobius transformations Kj. By one side we find 

1(1) -1(2) r l(3) r l(4) -1(5) r l(6)\ _ / r l(l) r l(2) r l(3) r l(4) -1(5) r l(6) 



-rr r 11 - 1 ' r 1 ^' r 11 -- 3 ' r 1 w r 1 ^ ' r-'-w — r 11 - 1 ' r 1 ^' r 1 ^' r 1 ^ r 11 - ' r 1 

A 13 I ^1 , "I !^1 !«1 ,V 1 I \ 1 ' 1 ' 1 ' 1 ' 1 ' 1 

that is, the transformations K\ (identity) and K13 of the first column of (|26p arc equivalent, a 
fact that holds for the transformations of any column. Thus we have to take into account only 
the twelve transformations of the first row. However, we find as well that 

k (r 1 ^ r 1 ^ r 1 ^ r 1 ^ r 1 ^) r^ 6 A - ( c 1 ^ r 1 ^ r 1 ^ r 1 ^ r 1 ^) r 1 ^ 

-ft 7 I «1 1^1 »V3l j^l I — I W ,^1 , ^1 ,^1 , ^1 I 
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that is, K 7 rearranges the initial kernels. This could be expected because K\ H > K 7 corresponds 
to the inversions (x, y) <-> (1/x, 1/y) which are already incorporated into the original set. This is 
the same for the other transformations corresponding to inverted mappings, that is, Ki •<=>• 
(i = 1, • • • ,6). Therefore, we may apply only the transformations K 2 , K 3 , K4, K§ and Kq on 
G^ according to Eq. ([25]) in order to find new sets, each one containing six hypergeometric 
functions. To determine which are the transformations Ni that are suitable to generate new 
kernels, it is sufficient to regard the pair ^G^ 1 ' , G^ 2 '^ ; the other pairs may be obtained by 

replacing the F (a, b; c;u) which appears in G 1 ^ by the hypergeometric functions written in 
Eqs. (|3TbT cV 

Thus, from K2 wc find 



Gl (1 \x : y) = \(l- X )(l-y)}- a F 



a, 1 + a — 5; 7; 



(a - l)xy 
a(x-l)(y-l) 



Gl {2) (x,y) = [(l-x)(l-y)} 



7 — O — 1 



(xy) 



1-7. 



(36) 



F 



together with the two pairs generated by using the hypergeometric functions (|31bl -c). as ex- 
plained above. For this case, the generators N2 and N3 are equivalent to each other and, 
consequently, the Ni afford only three additional sets. We find 

#iG? W o N 4 Gl {l \ iV 2 G 2W 7V 3 G 2W , iV 5 G 2(i) ^ JV 8 G£ (i) , AT 6 G! (i) N 7 G{ {l) . 

Thence wc may use only the transformations N2, N§ and Nq. On the other hand, by using K3, 
we get 



GfW) = 



f - 



F 



a, 7 + 5 - /3;7; 



(1 - a)xy 
(a - x)(a - y) 



G? (2) (x,?/) = [xy] 



1-7 



1-- 1-i 



7 — 1 — a 



(37) 



F 



a + 1- 7, 5+1-/3; 2- 7; 



(1 - a)xy 
(a - x)(a - y) 



This time iVi <^ 7V 6 , N 2 <^ iV 5 , N 8 and iV 4 iV 7 . Thus, the transformations N 2 , N 3 

and N4 are sufficient. The transformation K4 yields 



Gt W (x,y)=F 



a,0;6: 



(x-l)(y-l) 



I- a 



Gf 2 \x,y) = [{l-x){l-y)\ 



is 



F 



a+l-5,P + l-5;2-5: 



(x-l)(y-l) 



1 - a 



(38) 



Since Ni <^> N 3 , N 2 <^> N 4 , N 5 N 7 and N e <=> N 8 , we can choose only N 2 , N 5 and N e . By 
we get 





o - sx 


— a 




F 









a, 7 + <5 - /3; <5; 



a(x- 1)0,-1) 



Gf ) (x )2/ ) = [(l-x)(l- y )]^[(l-^) (l-|) 



a + 5, 7 + 1-/3; 2- 5: 



(x - a)(y - a) 



a(x-l)(y-l) 



(39) 



(x - a)(y - a) 



together with the sets obtained by applying N 2 , N 3 and N4, because N\ o- N 7 , N 2 O- N$, 
N 3 N$ and N4 Nq for this subgroup. Finally, the transformation K§ leads to 



G 6 ^\x,y)=F 



a,f3;e: 



(x - a)(y - a) 



a(a — 1) 



1-- 1-- 



F 



j + 5-a,^ + 5- f3;2-e; 



(x - a)(y - a) 
a(a — 1) 



(40) 
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and the sets resulting from this by A^, A3 and N4 since now N% A/5, N2 ATg, N3 •<=> N7 
and N 4 N 8 . 

In summary, the substitutions of the independent variables that convert the kernel equation 
© into a hypergeometric equation are (x,y) H> u(x,y), where u(x,y) is given by 



(l-a)xy (x-l)(y-l) a(x-l)(y-l) (x - a)(y - a) 



xy (a - l)xy 

a ' a(x — l)(y — 1) ' (a — x)(a — y) ' 1 — a (x — a)(y — a) ' a(a — 1) 

In general it is necessary to perform a substitution of the dependent variable as well. 



(41) 



2.4. Generalisation of Erdelyi's kernels 



These are given by products of hypergeometric functions containing an arbitrary separation 
constant A. When this constant is appropriately chosen, we recover the Lambe- Ward-type 
kernels. The Erdelyi kernels are obtained by rewriting equation ([5J in terms of the independent 
variables H 



a 



c 



(x - a)(y - a) 



(42) 



(1 - a)(xy - a)' 

and, then, by accomplishing separation of variables in the resulting equation. It turns out that 



G(s,y) = (!-£) A P(0 Q(C), 



(43) 



where P(£) and Q(C) satisfy the hypergeometric equation (Q} with the following sets of param- 
eters: 



P(£):a = a-A, b = j3 — X, c = 7; 

Q(C) :a = A, b = a + /3 — 7 — A, c = e. 

To show this, firstly we perform the substitutions (|42[) in Eq. (j4|. We get the equation 



(44) 



3^ G 

(i-oUa-o^ + 7-(a+/?+i)e 



C(i-C)0 + 



^ = 0, 

9C ' 



which, by the separation of variables G(£, £) = P(C)Q(C)i becomes 



(i-O 
p 



dP 



a/3P 



d 2 Q 



M C(1 - C) ^ 4 



0. 



Denoting the separation constant by A(a + j3 — 7 — A), we obtain 

dQ 



d( 



A(a + ^- 7 -A)Q = 0, 



(45) 



aP 



a/3 - 



A(a + (3 - 7 - A) 



P = 0. 



The additional substitution P = (1 — £) P leads to 



a" 2 P 



^ 2 



^(1 - Orirt + 7 -(a + (3 + 1- 2A)£ -7F - (« - - A)P = 



dP 



(46) 



In this manner, P(£) and <5(C) satisfy hypergeometric equations with the parameters given 
in (|4"4")) . If A = 0, we can take Q(C) = constant in order to recover the Lambe-Ward 
kernels, G} . On the other hand, for A = a and P(£) = constant, the resulting ker- 
nels belong to the generalised Lambe-Ward kernels G^\x,y) which accompany the kernels 
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(|39|) . since by relation (|33|) the arguments of the hypergeometric functions take the form 
u = (x — a)(y — a)/[a(x — l)(y — 1)], 1 — u or 1/u. 

For P(£) we select two hypergeometric functions in the neighbourhood of each singular point, 
as in the Lambe-Ward case. Since for Q(() there are six possibilities as well, we write the initial 
set as 



Gi^(x )tf ) = (l-f)- A pW(0 QW(C), 



(47) 



where pW(£) and Q (/) (C) 

are obtained from the six hypergeometric functions (|31al -c). having 
parameters specified in (|44[) . Explicitly, we find 



Ji (a ' ,) (x, tf ) = (*») 1 ^(l-f)" 



F 



a+l-7-A,/3 + l- 7-A;2-7; 



(48a) 



Q (0 (C); 



Gl (3 '°(x, 2/ ) = (i 



a / 



a-A,/3-A;l + a + ^- 7-2A;l 



'Hi 

a 



Q (i) (C), 



.1(4,0 



7 — a— /3+A 



F 



7-a + A,7-/3 + A;l + 7 - a- /3 + 2A;l-^l Q« (C); 

a 



(48b) 



Gi^(x,y) = (^-(l-f) 
Gi^W) = (^ (l-f) 



-A 



F 



F 



a — A, a + 1 — 7 — A:l + a — /?; — 

xy. 



/3-A,/3 + l-7-A;l + y S-a; 



Q (0 (C), 
Q (i) (C), 



where 

Q«(C) = F 



A, a + P - 7 - A; e; 



(x - a)(y - a) 



(1 - a)(an/ - a) 



(48c) 



Q (2) (C) = 



(x — a)(y — a) 



(1 - a)(xy - a) 



F 



(x - a)(y - a) 
(1 - a)(xy - a) 



Q (3) (C) = P 



A, a + j3 - 7 - A; 5; 1 



(x - a)(y - a) 
(1 - a)(xy - a) 



qW(C) 



(x - a)(y - a) 
(1 - a)(xy - a) 



1-5 



F 



e - A, 1 + A - 5: 2 - 5: 1 



(x - a)(y - a) 
(1 -a)(xt/-a) 



Q (5) (C) = 



(x — a)(y — a) 
(1 - a)(xy - a) 



1 -A 



A, 1 + A - e; 1 + 2A - a - /3: 



(1 - a)(xy - a) 
(x - a)(y - a) 



Q {6) (0 



(x - a)(y - a) 



(1 - a)(xy - a) 



A+7— a— 5 



x F 



(5-A,a + /3-7 + A;l-2A + a + /3: 



(1 - a)(xy - a) 



(x - a)(y - a) 



(48d) 



We can generate additional kernels with the same arguments for the hypergeometric functions 
by applying the index transformations iVj. However, as in the case of the Lambe-Ward kernels, 
we find that 

NiG{ {k,l) oN 2 G{ (k ' l \ AT 3 Gi (M) ^ NiG\ ik ' l) , N 5 G{ {k,l) -o- N 6 G\ (k ' l) , N 7 G\ [k ' l) & N 8 G\ ik ' l) . 
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Thence, again it is sufficient to use the transformations N$, N§ and N7 in order to generate 
the first subgroup of kernels. 

The next step refers to the generalisation of the Erdelyi kernels by means of the Mobius trans- 
formations Ki, which now change as well the arguments of the two hypergeometric functions. 
We transform only the kernel G\ whose explicit form is 



G^'foy) = (1- 



xy 



F 



a — A, /3 — A; 7 



xy 



x F 



A, a + (3 - 7 - A; e; 



(x - a)(y - a) 



(1 - a)(xy - a) 



(49) 



The 36 kernels are obtained by replacing each hypergeometric function by the other expressions 
given in Eqs. pial c). all of them with the same A. 

Again we can use only the transformations K2, K3, K±, K§ and Kq, the same ones employed 
in the Lambe-Ward case. In effect, -K13 gives 



K 13 G\ ilA \x,y) = (1- 



a J 



F 



a — A, P — A; 7 



x F 



A, a + /? — 7 — A; (5 



a(x-l)(y-l) 
(a - l)(xy - a) 



-1(1,3) 



(x,y)- 



Considering the other kernels, we conclude that K13 is equivalent to the identity K\ and, 
consequently, the kernels corresponding to the transformations of every column of matrix (|26|) 
are equivalent to each another. In addition, as in Lambe-Ward case, we find Ki 
(* = 1,--- ,6). 

Thus, applying K2, F3, K4, K§ and Kq on G* , we find the kernels 

— A 



G? (1 ' x) (x,y) = [(l-x)(l-y)]- a 



1 



(a - l)xy 



x F 



x F 



a — A, 1 + a — (5 — A; 7 
A, 1 + 2a - 7 - (5 - A; e; 

1 - 



o(x-l)(y-l). 
(a - l)xy 



a(x-l)(j/-l) 
(x - a)(y - a) 



-[(!-£)(!-«) 

x F a~\,^ + 5-/3 — A; 7 
x F X,a- j3 + 5 - X;S; 



a(l — x — y) + xy 

(1 - a)xy 
(x - a)(y- a) 

(1 - a)xy 
(a - x)(a - y) 
(x-l)(y-l) 
a + xy — x — y 



(50) 



(51) 



1 - 



(x-l)(y-l) 



1 - a 



-A 



F 



a — A, j3 — A; 7 



(x-l)(y-l) 



1 -a 



x F 



A, a — /3 — S — A; e; 



(x — a)(y — a) 



1 - 



a(a + xy ~ x — y) 



a(x-l)(y-l) 



(52) 



x F 
x F 



a — A, 7 + <5 — /? — A; <5: 



(x - a)(y - a) 
o(x-l)(y-l)" 



-A 



(x - o)(y - a) _ 



A, a — /3 + 7 — A; 7; 



xy — a 



(53) 



1 - 



(x - a)(y - a) 



1 -A 



a(a — 1) 



F 



a — A, /3 — A; e: 



(x - a)(y - a) 



a(a — 1) 



x F 



A, 7 + S - 1 - A; 5; 



a(x-l)(y-l) 
a(l — x — y) + xy 



(54) 
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The index transformations suitable for each set are: N 2 , N 5 and N 6 for G? (M) and G? (M) ; N 2 , 

N 3 and N 4 for G^' , Gf M) and G^ k ' l \ The previous kernels show that, in addition to the 
Erdclyi substitutions (|32")l for the independent variables, there are five other choices for (£,£)• 



3. Confluent Heun equation 

In this section we show that each transformation of the confluent Heun equation (CHE) is 
associated with a transformation of the equation for its kernels. We also show that some kernels 
of Heun's equation lead to kernels for the CHE through the limits <j5j) . For this end we rewrite 
the CHE © as 

d 2 d 
[M x — a] H(x) = 0, M x ^x(x-1)— +[--/ +{j + 8)x + px(x -!)]— + apx, (55) 

where in the previous equation M x is an ordinary differential operator. For the limits ([5]), the 
integral (J7J) becomes 

%{x) = e™ yi-\l-y) 5 - l G{x, y)H(y)dy, (56) 
-'in 

where the exponential results from the limit of (1 — y/a) e_1 . Besides this, Eq. ([8]) for G(x,y) 
takes the form 

[M x -M y ]G(x,y)=Q, (57) 
while the expression ([9]) for the bilinear concomitant becomes 

P(x,y) = e™ y~<(l-y) 5 



dG(x,y) ,dH(y) 
H (y> — 5 G(x,y) 



dy ' dy 



(58) 

First we discuss the transformations for the CHE (|5"5j) and its kernels (f5"T|) and, after this, we 
explain how to get kernels for the CHE from the ones of the general equation (JTJ) . 



3.1. Transformation of the confluent equation and its kernels 

There are 16 variables substitutions which preserve the form of the CHE Q. If H{x) = 
H(a] p, a, 7, 5; x) denotes one solution of the CHE (|55l) . these transformations are summarised 
in the rules T 1; T 2 , T 3 and T4 that operate as 

TxH(x) = (l-x) 1 ' 5 H[a-j{l -6);p,a+ 1-5,7,2 -6;x], 
T 2 H(x) =x 1 -~<H[<j+ (I-7) (p-S);p,a + 1 - 7, 2 -j,S;x], 

(59) 

T 3 H(x) = e px H[a - jp; -p,-y + S- a, 7, 5; x], 
T4H (x) = H[a — pa; —p, a, S, 7; 1 — x] . 

Compositions of these give the group having 16 elements. 

On the other side, equation (|57|) for the kernels is written in terms of the same differ- 
ential operator which appears in the CHE ([6]). Consequently, proceeding as in the case 
of the general Heun equation, the corresponding rules Ki for transforming a given kernel 
G(x,y) = G(p,a,7,6;x,y)) are 

K 1 G(x,y)= [(1 - x)(l - y)] 1 - 5 G [p, a + 1 - 5, 7, 2 - 5; x, y] , 
K 2 G(x, y) = (xy) 1 " 7 G [p, a + 1 - 7, 2 - 7, 6; x, y] , 

(60) 

K 3 G(x,y)=e-p^ G[-p, 7 + 6 - a, 7, 6; x,y], 
K 4 G(x, y) = G [ - p, a, 5, 7; 1 - x, 1 - y] . 

These rules can be verified by substitutions of variables. However, they are useful to produce 
new kernels when one knows an initial kernel for the CHE. In the following we obtain initial 
kernels as limits of kernels for the Heun equation. Any kernel which can be generated by the 
above transformations is omitted. 
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3.2. "Lambe- Ward-type" kernels 

The limits (J5)) applied to kernels of Sec. 2.3 give three kinds of kernels for the CHE, two 
in terms of confluent hypergeomctric functions and one in terms of Gauss hypergcomctric 
functions. The regular and irregular confluent hypergcomctric functions, denoted by $(a, c;y) 
and ^(a, c;y), respectively, are solutions of the equation 

d 2 ip , «, dip . „ . 

^ + (C - u) £- a ^°- (61a) 

The following types of solutions for Eq. (|61al) 

ipW(u) = $(a, c; u), ipW(u) = e u u 1 - c <S>{l - a, 2 - c; -u), 
ip( 3 )(u) = *( a ,c;u), <pW(u) = e u u 1 ~ c ^(l -a,2-c; -u), 



(61b) 



are all defined and distinct only if c is not an integer [l4[ ■ Alternative forms for these solutions 
follow from the relations 

$(a, c;u) = e"<f>(c - a, c; -u), *(a, c; u) = u 1 ^ c ^(l + a - c, 2 - c; it). (61c) 

In the present context, the above confluent hypergcomctric functions result from the limits 



1. 



lim F ( a,b; c; 1 — — ") = lim F (a, b; c; — — ") = i^Vt^a, a + 1 - b; u), 



US c— >oo 



lim F 

b— ^C30 



(a,b;c; ^ = $(a,c;«). 



(62) 



Sometimes, before applying these limits, it is necessary to use the relations (|3"2"T) and/or (|33[) 
and, in addition, multiply the kernels by suitable constants depending on the parameter a. 
Furthermore, note that both a and /? tend to infinity but such that (3 = —pa, where p is 
constant. For this reason we can write, for example, 



F 



pxy 

/3 



and, thence, accomplish the limit j3 — > oo by keeping pxy hxed. 
Thus, the first set of section 2.3, Gj (x,y), gives the limits 

Gj (1) (i,!/)^$ (a, 7; -pxy) , 
G^ (2) (x, y) -> e-*"* (xy) 1 " 7 $ (1 - a, 2 - 7; pxy) , 
Gi (3) (z, y) and G 1 (5) (x, y) -> * (a, 7; -pxy) , 
Gj (4) (x,y) and Gj (6) (x,y) -> e-^(xy) 1 -T * (1 - a, 2 - 7; pxy) , 
which, by means of (|61b[) . can be written as 

G < ( ) (x,y) = ip { - i) {u), with u = -pxy 7 a = a, c = 7 [i = 1,2,3,4]. (63) 

This initial set of kernels is really due to Lambe and Ward because it is given by confluent 
hypergeometric functions whose arguments depend on the product xy as in [[J. The previous 
transformations produce new solutions of this type. 

The limits of the set G^ (x, y), Sec. 2.3, gives six kernels in terms of Gauss hypergeometric 
functions (|31al -c). written as 

G^(x,y) = [(l-x)(l-y)]- a F^(u), [* = 1, • - ■ ,6] (64) 

with u — r-r -, a = a, b = 1 + a — 8, a = 7. 

(x- l)(y- 1) 

We can also apply the transformations Kj in order to generate a group of kernels. Notice that 
this group results from a generalisation of the Lambe- Ward kernels and, as far as we know, it 
is new. 
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Finally, another set of kernels is given by confluent hypcrgeometric functions f|61b[) whose 
arguments are proportional to x + y — 1 . It is obtained as limits of the kernels f)40[) and reads 

&f =ip { - i) {u), withu = -p(x + y-l), a = a, c = 7 + 5 [i = 1,2,3,4]. (65) 

Other kernels follow from the transformations Ki . This group is also a result of a generalisation 
of the Lambe-Ward kernels, but kernels having this form are already known in the literature 

3.3. "Erdelyi-type" kernels 

By taking the limits of the Erdelyi-type kernels we find two groups of kernels for the CHE. 
The group given by products of confluent hypergeometric functions has already appeared in 
the literature [l5[, whereas the group given by products of hypergeometric and confluent hy- 
pergeometric functions seems to be new. In the limit process we suppose that A is kept fixed, 
that is, we assume that A does not depend on the parameters /3 and a of the Heun equation. 

Thus, by taking the limit of the kernels G 1 ( " fc ' Z ' ) (x, y) given in (|47[) . we get an initial set of 
kernels given by 

Gf J) =^)(0^')(0, [i,j = 1,2,3,4] (66) 

where <pW (£) and (pffl (£) are the solutions (|61b[) for the confluent hypcrgeometric equation, 
having the following arguments and parameters : 

<p l (0 : £ = -pxy, a = a - A, c = 7; tp J (() : C = pip - l)(y - 1), a = A, c = 5. (67) 

The four kernels given by products of regular functions $ are 

G< M) (x, y) = $ [a - A, 7; -pxy] $ [A, 6; p(x - l)(y - 1)] , 

G?'*\x,y) = e^-^y^[{x-l){y^l)] l - 5 ^[a-\, T ,-pxy] 

x $[l-\,2-5]-p(x-l)(y-l)], 
G^ l \x, y) = e-P^ixy) 1 -^ [A + 1 - a, 2 - 7; pxy] $ [A, 6; p(x - l)(y - 1)] , 
G?' 2 \x,y) = e^ + ^(^) 1 ^[(x-l)( 2 /-l)] 1 - 5 <I>[A + l-a,2-7;p^] 



x $[l-\,2-5;-p{x-l)(y-l)}. 

Replacing one or both <!> by \& we get the set with 16 kernels. The other sets result from this 
by the transformations Kj. 

In the second place, the kernel G^ 1 ' 1 ^ (x, y), Eq. (|54|) . yields a kernel G < f ,1 \x^y) constituted 
by a product of hypcrgeometric and confluent hypcrgeometric functions, namely, 

G[ 2A) (x,y) = (l-x-y)- X ^[a~X, 1 + 5-2X;p(l-x-y)} 



x F 



A, 7 + S - 1 - A; 5: 



(x-l)(y-l) 



l-x-y 

By considering the limits of the full set G^ (x, y) we obtain the initial set 
(x, y) = <p® (OF^ (C), [i = 1, . . . , 4; 3 = 1, • ■ • , 6] 



(68) 



where ip^ are the four solutions for confluent hypergeometric equation and where Fv) are the 
six solutions for hypergeometric equation with the following arguments and parameters: 

£ = P(l-x-y), a = a-A, c = 7 + (5-2A 

(69) 

F U) (C)-- C= ^-i-^ . a = A > b = 7 + 5-l-A, c = <5. 



The kernels G^ (a;,y), Eq. (j5T))) . also lead to kernels given by products hypergeometric and 
confluent hypcrgeometric functions, but we may show that these are connected with the previous 
ones by the transformation K4. 
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4. Concluding remarks 

We have taken the following steps to deal with kernels for integral relations among solutions 
of Hcun equations: 

• the use of an integral with a weight function w(x, y) which allows to write a given Hcun 
equation and the respective equation for the kernels in terms of operators functionally 
identical, say, M x and M y ; 

• the use of the known transformations of the equation in order to get the actual form for 
the transformations of the kernel equation; 

• the generation of new kernels by applying the previous transformations to an initial kernel 
or set of kernels. 

• the use of a limiting procedure to generate kernels for the confluent Hcun equation. 

For the (general) Heun equation we have used Maier's transformations, discarding as inappro- 
priate the forms given in Refs. [f| and [l2j for the homotopic transformations. As initial kernels 
we have employed the ones found by Lambe and Ward, Eq. and by Erdelyi, Eq. (|4"T)) . 

In this manner, in section 2 the transformations for integral relations have afforded several 
new kernels for the Heun equation, given by a single hypergeometric function and by products of 
two hypergeometric functions. We have seen that only six of the homographic transformations 
for the kernels are effective, namely: K±, K 2l - ■ • ,Kq, where K\ is the identity. The fact these 
are just the six first transformations of Appendix B is a consequence of manner in which we 
have written the elements of matrix (fTU)) . 

We have written only some of the possible kernels, but a wider list can be generated by 
index transformations which lead to new kernels where the hypergeometric functions possess 
the same argument but different parameters. In addition, from a kernel with a given argument, 
new kernels follow from the fact that the hypergeometric equation formally admits solutions 
with different arguments in the vicinity of each singular point. 

In section 3, the confluence procedure (O has led to five sets of initial kernels for the CHE, 
three of them arising from generalisations of the Lambe- Ward and Erdelyi kernels by means 
of Mobius transformations. This in association with the Leaver version for the CHE and the 
concept of Whittaker-Ince limit suggest new kernels also for the double-confluent Heun equation 
(DCHE) and for limiting cases of the CHE and DCHE. In effect, by substitutions of variables, 
the CHE © can be written in the Leaver form namely, 

z(z-z )^ + (B 1 +B 2 z)^+[B 3 -2r ] u}(z-Zo)+u> 2 z(z-z )]U = 0, [u + 0] (70) 

where Bi, rj and uj are constants, and z = and z = zq are the regular singular points. When 
zq = (Leaver's limit), this gives the DCHE 

z2 ^ + ( B i+ B 2 z )^ + (, B 3- 2r l L < J z + u 2 z 2 )U = 0, [Bi^O, u^O], (71) 

where now z = and z = oo are irregular singularities. On the other side, these equations 
admit the limit [lj| 

uj — > 0, rj — > oo, such that 2i]ui = —q, [Whittaker-Incc limit] 

where q should not be confused with the parameter q of the Hcun equation ([T]). The Whittakcr- 
Ince limit of the CHE and DCHE are, respectively, the equations 

z(z-Zo)^ + (B 1 +B 2 z)^ + [B 3 +q(z-z )}U = 0, [q^O] (72) 

~dz^ + (Sl + B2Z) ~dz~ + + qZ) U = ° ! [q * °' Bl * ° ] (73) 

which have a different type of singularity at x = oo as compared with the original CHE and 
DCHE [l6l - [l8| . As the preceding limits connect the CHE with the DCHE and their respective 
Whittaker-Ince limits, we may expect to find kernels for each of these out of kernels arising from 
the Hcun equation. Thus, by developing the results of section 3 we could unify the treatment of 
these equations. We advance that we will find that the usual kernels for the Mathieu equation 
(20j are particular cases of the kernels for Eq. (|T2")) . 

The construction of new kernels can be envisaged as a first step for seeking new solutions for 
the Heun equations by means of integral relations [l9j]. However, as in the case of the Mathieu 
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equation [20J, it is not easy to use this technique. Thus, in the following we focus on possibility 
of using some of the Maier transformations to extend the solutions in series of hypergeometric 
functions given by Svartholm in 1939 [2l| and by Erdclyi in 1944 [l3j]. Also in this context, the 
previous limits become relevant. 

As a prior consideration, we note that there are three types of recurrence relations for the 
series coefficients - not just one as given in the original articles and repeated since then [H, 12Tl. 
The two additional relations may be found by the procedure used in Appendix A of Ref. 22[ . 
Then, by virtue of the new relations and by means of the homotopic transformations, we may 
show that the Svartholm solutions include as particular cases the eight Fourier-type solutions 
found by Ince in 1940 for the Lame equation [l2l [23j . 

On the other side, both the Svartholm and the Erdelyi solutions are valid only if the parameter 
a satisfy the condition Re-^/1 — (1/a) > which is assured by requiring that a ^ [0,1]. In 
addition, the former solution is given by a single expansion in series of hypergeometric functions 
and converges only over a finite region of the complex plane; however, the latter is in fact a 
set of expansions in terms of hypergeometric functions which, by analytical continuation, may 
cover the entire complex plane provided that a characteristic equation is fulfilled. 

Then, Erdelyi's solutions are candidates to solve a cosmological problem formulated by Kan- 
towski |24| because in this case the variable x extend to infinity. Nevertheless, the problem 
demands solutions valid also for a 6 [0, 1]. These may be derived by applying on the Erdclyi 
solutions one of the following linear transformations: 

M 9 H(x) = H Q, 0, 7 , e; f ) , M 61 H(x) = H ^s£; a, 0, 6, e; ff±) , 
M 101 H(x) =H(^, ^±^;a, 0, e, 7 ; , 
M lm H(x) =H(^ T , ^±M£; a, 0, e, 6; 

Hence we find, respectively, the conditions: a [1, oo), a ^ (— oo, 0], a 0[1, oo) and a $ (-co, 0]. 
However, for each case it is necessary to reexamine the domains of convergence. 

After these preliminaries, we conclude by adding that some of the Svartholm and Erdelyi 
solutions - the solution (|21[) , for example - lead to solutions for the CHE by means of the limits 
((5|). To prove this it is sufficient to divide the recurrence relations by the parameter a before 
performing the limits. Furthermore, we can show as well that the solutions for the CHE in 
the form ([70]) supply solutions for the DCHE (fTTj) through the Leaver limit (zq — > 0). These 
are additional reasons for choosing the Svartholm and Erdelyi solutions as a starting point for 
further investigation. 



Appendix A. Equations of the first section 

Equations ((HJ and (j9|) as well as the condition P(x,yi) = P(x,y2) are obtained from the 
general theory of integral relations which is established for w(x,y) = 1, that is, for 

K(x) = r K(x,y)H(y)dy (Al) 
Jyi 

where K(x, y) denotes the kernel. In this case the equation for K(x, y) is given in terms of the 
operators M x and M v , where M y is the adjoint operator Q corresponding to M y , that is, 

M y = y(y - l)(y - a)— + [(2 - j)(y - l)(y - a) + (2 - 8)y(y - a) + (2 - e)y(y - 1)] — 
+ [4 - 2(a + + 1) + a/3] y + 0(7 + 6 - 2) + e + 7 - 2. 



By applying M x to the integral (|Al|) 

and supposing that the integration endpoints are independent of x, we find 



hi 



M x X(x) = / H(y) [M x - M y ] K(x, y)dy + / H(y)M y K(x, y)dy. 
Jyi Jyi 

Then, by requiring that the kernel satisfies the partial differential equation 

[M x - M y ] K(x, y) = 0, (A2) 
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the right side of of the previous integral reduces to J H(y)M y K(x,y)dy. Using the Lagrange 
identity 

H(y)M y K(x, y) = K(x, y)M y H(y) + |-P(x, y) © qK(x, y)H{y) + |-P(x, y), 



H(y)^l-K( X ,y) dH ^ 



where now P(x, y) is given by 

P(x,y) = y{y - l)(y - a) . 

ay ay 
+ [(1 - 1){V ~ 1)(» -a) + (l- 5)y(y - a) + (1 - e 

then we find that the integral yields 



M x 3i(x) 



.72 



Hi 



qK(x,y)H(y) 



dP(x,y) 
dy 



dy&[M x -qMx) = P(x,y) 



(A4) 



Therefore, %{x) will be a solution of the Heun equation if K is solution of (|A2j) . if the integral 
(|A1|) exists and the limits of integration are so chosen that P(x,yi) = P(x,y2)- Further, by 
setting WL(x,y) = w(x,y)G(x,y), where w(x,y) is defined in Eq. ([7]), we recover equations ([5]) 
and ([9]). Notice that Eq. (|A2[) is inadequate to deal with the kernel transformations because 
the operators M x and M y do not present the same functional form. 



Appendix B. Mobius transformations for kernels 

As in the index transformations, firstly we find the Mobius transformations Mi for H(x) in 
Maier's table and, then, write the kernel transformation Ki in accordance with the rule (f2"4")) . 
Thus, the 24 expressions for the Ki of matrix (f2"5)) are the ones given below. 



KiG{x,y) = G(x,y) = G [a; a, 0, 7, 5; x,y] , [Identity]. 



K 2 G(x,y) = [(l-x)(l-y)]- a G 



a x y 

' a, 1 + a — 0, 7, 1 + a — p; ■ 



a - 1 



x - 1 y - 1 



•^3G(x, y) = 



1-^ 

a/ V a 



■ a, 7 + 7,1 + a 



1 -a 



a; — a y — a 



K 4 G(x, y) = G [1 - a; a, /3, <5, 7; 1 - x, 1 - y] 



^5G(x, y) = 



1-^ 

aJ \ a 



1 x — 1 y — 1 

-; a, 7 + 6 - 0, 6, 1 + a - 0; , 

a x — a y — a 



K e G(x, y) = G 



a— 1 a — x a — y 
;a,0,e,r, , 



K 7 G(x,y) = (xy) a G 



-; a, 1 + a — 7, 1 + a — j8, 5; — , - 
a x y 



K 8 G(x,y) = (xy) a G 



a ~ 1 -, r-, n ^ - 1 2/ ~ 1 
; a, 1 + a - 7, 0, 1 + a — 0; , 

a x y 



K 9 G(x,y) = (xy) a G 



1 — a; a, 1 + a — 7, e, 1 + a — 0; 



x — a y — a 
x ' y 



K 10 G(x,y) = [(l-x)(l-y)}- a G 



a, 1 + a — 5, 1 + a — 0, 7; ■ 



1 -a 



1 — x 1 — y 
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K 11 G(x 1 y) = [(l-x)(l-y)]- a G 
K 12 G(x,y) = 
Ki 3 G(x,y) = G 



a; a, 1 + a — 5, e, 1 + a — /3 



x — a y — a 



x — 1 ' y — 1 







— a 




G 









a a a 

a, 7 + o — p, 1 + a — p, 7; 



a - 1 



a — x a — y 



1 a x V 

-;u,p,7,e, -, - 

a a a 



K u G(x,y) = [(l-x)(l-y)]- a G 



K 15 G(x,y) 



K 16 G(x,y) = G 



K 17 G{x,y) 



K 1& G{x,y) = G 



1 - 



a - 1 (a-l)x (a - l)y 

; a, 1 + a - (5, 7, e; — — , — — 

a a(x — 1) a{y — 1) 

(1 - a)x (1 - a)y 



1 — a;a,7 + <5 — /?,7,5; 



2; — a y — a 



1 — a 1 — a 



1 -a 

x 
a 



a; a, 7 + 8 — p\ J, 7 



a{x - 1) a(y - 1) 
x — a ' y — a 



a a a-x a-y 



a - 1 



ATi 9 G(x,2/) = (xy) " G 



A" 20 G(x,?/) = (xy) Q G 



a — 1 a — 1 



a; a, 1 + a — 7, 1 + a — (3, e; 



a a 



a a(x-l) a(y-l) 

a, 1 + a — 7, 0, e; 



a- 1 



(a - l)x' (a - l)y 



1 -a 



2f 2 iG(x,y) = (xy) a G 
K 22 G(x,y) = [(l-x)(l-y)r a G 
K 23 G(x,y) = [(l-x)(l-y)]- a G 
K 24 G(x, y) 



1 x-a y-a 

■ a, 1 + a — 7, e, 0; ■ 



(1 - a)x' (1 - a)y 



1 — a;a,l + a — (5, 1 + Q 



1 — a 1 — a 



1-x 1-y 







— a 


[('-;) 


G 









1 x-a y-a 

-; a, 1 + a - d, e, 7; — — , — — 

a a(x — lj a(y — 1) 



a - I r. „ r 1 — a 1 — a 
; a, 7 + — p, 1 + a — p, 



x — a y — a 
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